Failure of geometric electromagnetism in the adiabatic vector Kepler problem 
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The magnetic moment of a particle orbiting a straight current-carrying wire may precess rapidly 
enough in the wire's magnetic field to justify an adiabatic approximation, eliminating the rapid 
time dependence of the magnetic moment and leaving only the particle position as a slow degree 
of freedom. To zeroth order in the adiabatic expansion, the orbits of the particle in the plane 
perpendicular to the wire are Keplerian ellipses. Higher order post-adiabatic corrections make the 
orbits precess, but recent analysis of this 'vector Kepler problem' has shown that the effective 
Hamiltonian incorporating a post-adiabatic scalar potential ('geometric electromagnetism') fails to 
predict the precession correctly, while a heuristic alternative succeeds. In this paper we resolve the 
apparent failure of the post-adiabatic approximation, by pointing out that the correct second-order 
analysis produces a third Hamiltonian, in which geometric electromagnetism is supplemented by a 
tensor potential. The heuristic Hamiltonian of Schmiedmayer and Scrinzi is then shown to be a 
canonical transformation of the correct adiabatic Hamiltonian, to second order. The transformation 
has the important advantage of removing a 1/r 3 singularity which is an artifact of the adiabatic 
approximation. 



I: INTRODUCTION 

An analogue to the Kepler problem of motion in an 
inverse square force can be realized with magnetism and 
adiabaticity - except that the analogue to the planetary 
mass is not a scalar, but a vector component jy. For a 
particle of mass M with spin s and gyromagnetic ratio 
/i, moving in a magnetic field B{x), the Hamiltonian is 



2M P 



(1) 



If the field B = &9r~ x is due to a uniform current flowing 
in a straight, thin wire along the z-axis, the z component 
of linear momentum is conserved, and so the problem 
reduces to finding the orbit of the particle in the xy plane. 
The magnitude s of s is fixed, and the z component of 
total angular momentum J z = xp y — yp x +s z is a constant 
of the motion. In realistic cases the dimensionless ratio 
e = sJJ z is small. This implies a separation of time 
scales, so that adiabatic methods may be applied pi 0, 

To zeroth order in e one can (as we will review below) 
replace H with the effective Hamiltonian 
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where we drop the z direction, rescale to dimensionless 
units, and use standard polar co-ordinates r, 9, centred 
on the current-carrying wire. In Hado the pairs r, p r 
and 0,pq are canonically conjugate, and a is simply a 
constant. We can recognize Hado as the Hamiltonian 
for Kepler's problem of motion under an inverse square 
force. See Figure 1. In this vector version of the prob- 
lem, however, the analogue of the particle mass, a cx se, 
is a vector component. It may therefore be positive or 
negative; but since the problems considered in this paper 



are much less salient for unbound motion, we will assume 
a > 0. 
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FIG. 1: Orbits of the magnetic particle about the current- 
carrying wire at the origin, from numerical solution of the 
equations of motion (|To[) - (|ll|) derived from (|l|). Initial con- 
ditions are a(0) = 0.9, a(0) = 0, r(0) = 10, p r (0) = 0, and 
e — s/Jz — 0.2. The slightly polygonal appearance of the 
trajectory, with rounded corners between straighter sections, 
is not a numerical artifact: it is the rapid component of the 
particle motion. 

As Figure 1 illustrates, for small e the particle orbits 
are indeed well described as ellipses. Two kinds of correc- 
tions, however, are also evident: a rapid wobble, which is 
straightforward to compute, and a slow precession, which 
is not. Capturing effects like this slow precession is the 
goal of higher order adiabatic approximations, which add 
to Hado terms of higher order in e. Classic papers on 
adiabatic methods ||, ||[ ||, ^, |^ prescribe for this problem 
the second order effective Hamiltonian 



H, 



GEM 
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whose correction term can (as we will explain) be called 
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'geometric electromagnetism'. But see Figures 2 and 3: 
Hgem is actually worse than Hado, in that it yields 
precession in the wrong direction; whereas the precession 
is given with impressive accuracy by a slightly different 
alternative that has been identified, without derivation, 
by Schmiedmayer and Scrinzi M: 
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So what has gone wrong with adiabatic theory? And 
what is Hsst 
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FIG. 2: Fhe same as Fig. 1 (solid curve), with the geometric 
electromagnetism approximation superimposed for compari- 
son (dashed curve). The dashed curve is the numerical so- 
lution of the canonical equations of motion derived from the 
Hamiltonian Hgem of (H), with the same initial conditions 
as in Fig. 1. 
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FIG. 3: Again the same as Fig. 1 (solid curve), with the 
heuristic approximation of Schmiedmayer and Scrinzi super- 
imposed for comparison (dashed curve). The dashed curve is 
the numerical solution of the canonical equations of motion 
derived from the Hamiltonian Hss of (Q), with the same e 
and initial conditions as in Figs. 1 and 2. 

In this paper we explain that nothing has gone wrong 
with adiabatic theory. The problem is simply that a 
tensor potential, formally identified by Littlejohn and 
Weigert pf, and proportional in our case to ?'~ 3 , must 



be included in addition to the more widely known vec- 
tor and scalar potentials of geometric electromagnetism. 
In the next four Sections (II- V) we build up to a simple 
derivation of this term, first reviewing the zeroth order 
approximation, and then geometric electromagnetism. In 
Section IV we discuss the subtlety of relating the ini- 
tial conditions for the exact and approximate evolutions. 
Then in Section V we present the tensor potential for the 
classical version of our vector Kepler problem, and also 
provide a simple physical explanation of the phenomenon 
(which is much more general). In Section VI we show 
that Hss is the result of a canonical transformation of 
the correct second order adiabatic Hamiltonian; and we 
show that this transformation is particularly convenient 
because it removes the r -3 singularity. After treating 
the quantum version of the problem in Section VII, we 
discuss our results and conclude. 



II: THE VECTOR KEPLER PROBLEM VIA 
ADIABATIC APPROXIMATION 

Dropping the trivial p 2 term, the Hamiltonian (Q) may 
be written 



H = 



Pr 

2M 



Li 



B[L 



r 



9) (5) 



where re 
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2Mr 2 

iy defines the usual polar co-ordinates 
for the particle's position, p r = p x cos9 + p y sin9^ L z = 
xp y — yp x - We introduce the angle <fi which is canonically 
such that s T 



conjugate to s z 



V^ 2 



and 9, and p r and r, are of course the other canonically 
conjugate pairs. 

It is convenient to use polar spin axes, so that we have 
sg,s r instead of s x ,s y , and the canonical pair sg,a in- 
stead of s z , 4>. Defining the new canonical momentum 



SB 



s r sin i 



\sui(<f>-6) (6) 



mixes the angular co-ordinate 8 with the spin sector, 
though. Checking the Poisson brackets, we find that to 
keep the transformation canonical, we must also change 
the momentum conjugate to 9 from L z to J z . So we must 
write 



L z = J z 



= J Z - 



(7) 



It is also convenient to make all our variables dimen- 
sionless, by rescaling them in terms of Hamiltonian co- 
efficients and a reference angular momentum J: 

J 3 . „ M(sBii) 2 



t 



M(sBn) 2 
J 2 



t 
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J z — > Jpe 
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These are the variables with which we will proceed. Of 
course pg is a constant of the motion, and so we can (and 
do) always set the numerical value of pg — > 1 by choosing 
J = J z ; but we retain it as a canonical variable, which 
is needed in the equation of motion for 6. In terms of 
these variables the exact Hamiltonian (||) can be finally 
re- written 



H = 



P7 



(pg — eyl — a 2 cos a) 2 



a 
r 



(9) 



2 2r 2 

where e = s/ J. 

As we have defined it, a is of order e° in the cases we 
consider, which is convenient for keeping track of powers 
of e; but the dimensionless momentum canonically conju- 
gate to a would be sg/ J = ea. This gives the equations 
of motion for a and a an extra factor of e _1 : 



ldH r ~ VT 

= — (pg — e v 1 — & cos a) 
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Our other equations of motion are the usual canonical 
ones: 

. _ dH 

dp r 

dH (p g — eyi — a 2 cos a) 2 
dr r 3 
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We can now obtain the vector Kepler Hamiltonian (g) 
from (0) simply by setting e — > 0. In this limit, a evolves 



infinitely fast according to (10); but a plays no role in the 



Kepler analogy, and can simply be ignored. On the other 
hand, because a evolves so rapidly, from the equation for 
a in ( |ll|) we can easily see that 



a(t) 



spe- 



cos a 



0(e 2 ) 
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where a is a constant. So to zeroth order in e, H is indeed 
equivalent to Hado, a s claimed. And as Fig. 1 confirms, 
the evolution of (r, 6) under Hado does provide a close 
approximation to that under H, over intermediate time 
scales. This is not bad; but the errors in using Hado do 
become large for t of order e~ 2 . We can do better than 
this. 



Ill: POST-ADIABATIC CORRECTIONS: 
GEOMETRIC ELECTROMAGNETISM 

To improve on our simple adiabatic approximation we 
must become systematic about precisely how we are ex- 
panding in powers of e, because there is more than simple 



perturbation theory going on. The very rapid evolution 
of a is a warning of this. In fact the approximate con- 
stancy ('adiabatic invariance') of a, on which the Kepler 
analogy is based, is directly connected to the rapidity of 
a. The fact that a is fast means that the Hamiltonian 
varies steeply as a function of its conjugate momentum 
a, and this in turn means that, since H is identically con- 
served under time evolution, a must vary only slowly. So 
the rapid variation of a is no accident. And the essential 
feature of our adiabatic approximation is that as well as 
dropping terms whose magnitude is higher order in e, we 
are also dropping terms, like a, whose variation is on the 
rapid timescale e~ x t. 

An elegant and general treatment of post-adiabatic 
corrections is available in Ref. Q|. The result that is 
reached is physically intuitive: to eliminate rapid evo- 
lution while preserving slow evolution, one must remove 
rapid components from the Hamiltonian, by averaging 
over the rapid angular variable which is canonically con- 
jugate to the adiabatic invariant. If we identify a as the 
adiabatic invariant, then a is the rapid variable, and so 
we should integrate H over a. In many problems of this 
sort, it is common that rapid components of the Hamil- 
tonian appear, as they do in ([}]), as shifts in the posi- 
tional momentum. Thus to first order in the adiabatic 
small parameter (in our case, e), after integrating out 
the rapid co-ordinate one obtains an effective Hamilto- 
nian term linear in spatial momentum, A ■ p for some 
function A. In our case, this first order A vanishes, be- 
cause the integral over 2tt of cos a is zero. But in general 
such terms do not vanish. They imply that a neutral 
particle with a magnetic moment, moving in a spatially 
varying magnetic field, experiences a force as though it 
were a charged particle moving in an effective magnetic 
field, whose vector potential is A. This phenomenon has 
therefore been named 'geometric magnetism'. 

At second post-adiabatic order, one typically encoun- 
ters a scalar potential; in our case, integrating (||) over a 
produces the scalar term 



1 

27 



da 



e 2 (l-a 2 ) 



a 2 ) cos 2 a 
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which produces Hqem as given by (||). Continuing the 
analogy begun at first order, this post-adiabatic scalar 
potential may be compared to an effective electric poten- 
tial. The effective Hamiltonian produced by averaging H 
over a can therefore be said to include geometric electro- 
magnetism. 

So why does it not work very well? Since integrating 
out a has been quite straightforward, we might wonder 
whether the observed inaccuracy of Hqem is due simply 
to not going on to third or higher order in e. After all, 
e = 0.2 as in Figs. 1-3 is not really such a small parame- 
ter. But H contains no terms beyond order e 2 , so if the 
rule really is to average over a, then it is not clear where 
higher order corrections could come from. And in a case 
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where e is truly small, so that different orders in s are 
quantitatively obvious, we can check (see Figs. 4 and 5) 
that the discrepancy between H and Hqem is indeed a 
problem of lower order in e. In fact there are two inde- 
pendent problems, which we will discuss in turn, starting 
with the larger but more trivial one. 

e" 1 Ar 
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FIG. 4: The error Ar(t) = f-r and A0(t) = 9-9, where f , 9 
are r, 9 as given by the geometric magnetism approximation. 
The quantities plotted are obtained by numerically solving 
the same two sets of equations as in Fig. 2, but with e = 0.01 
and f(0) = r(0) = 1.5 (and all other parameters the same as 
in the previous Figures). The errors in the pure geometric 
electromagnetism approximation are of order e, rather than 
e 3 as desired. For these parameters, however, the heuristic 
approximation of Schmiedmayer and Scrinzi does not fare any 
better: its errors are not distinguishable to the eye from those 
shown here. 



IV: INITIAL CONDITIONS 

It is surprising that Fig. 4 shows errors of order e in 
evolution under both Hqem and Hss, since they agree 
with Had2 to that order; we should expect them to fail at 
0(e 2 ). In fact this expectation is correct, and the reason 
for the even poorer performance shown in Fig. 4 really 
has nothing to do with the form of the Hamiltonian. The 
problem is with using the same initial conditions for the 
adiabatically approximate evolution as for the exact evo- 
lution. 

As we can see in Eqn. (|l2|), the system's several vari- 
ables generally evolve as a superposition of slow and 
fast components (where the fast components may also 
be slowly modulated). Adiabatic approximations seek to 
reproduce only the slow components. This means that 
initial conditions need some consideration: the exact ini- 
tial values Pr(0), cr(0), etc., are the instantaneous sums 
of slow and fast components of motion; but the adia- 
batic approximation requires the initial values of the slow 
components alone. We must therefore identify the fast 



components in our variables, and subtract their instan- 
taneous values at t = from the exact initial conditions, 
to obtain the initial conditions that should be used with 
adiabatically approximate evolution. 

Fortunately it is not at all difficult to identify the fast 
components. If in ( pd| ) we examine the equation of mo- 
tion for 6, for instance, we see 



pe \A - c 2 
— — £ cos a ^ 



(14) 



Comparing this with the very rapid evolution of a ac- 
cording to (EOT) , it is obvious that we must have 



Bit) =8(t)+e 2 sma ^ /1 g -+0(e 3 ), (15) 
r 

where 8 is the slow component, which does not involve 
a. Similarly we find, up to corrections of order e , 



r(t) - f{t) 

p r (t) = p r (t) + E 2 



^/T^ 2 



■ sin a 



a(t) = a - ep g 



Vl — o 2 cos a e 2 apg 
r 2r 2 



(16) 



pePrV^ — v 2 sin a (1 — a 2 ) cos 2a 



4r 



where a is still exactly constant, and not just slow. Ac- 
tually, the e 2 apg/r 2 term in ait) is not trivial to deter- 
mine; but it can be obtained fairly easily using some of 
the insights of the next Section, and we include it here 
for completeness. 

So when we evolve the slow components under Had2, 
therefore, we should use r(0) — > f(0), p r (0) — > p r {0), 
etc., and a — > a, where the barred quantities are given 
in terms of the exact initial conditions by 



f(0) = r(0) 

Pr(0) = Pr{0)~e 2 



r 2 (0) 



sina(0) 



6»(0) = 6>(0) -e 2 sin a(0) 



y/l-^ 2 (0) 
r(0) 



yi-^Wcos^O) 2 a(0)p 2 (0) 
a = (7(1) j + e 7^ e 



"(0) 



2r 2 (0) 



; PePr(0)Vl-g 2 (0)sma(0) 
r(0) 

2 [1 - cr 2 (0)] cos 2a(0) 
4r(Q) ' 



(17) 



We might be puzzled at this point by the fact that in 
replacing c(0) with a we are actually throwing away some 
components that are non-constant but slow, rather than 
fast; this does not really seem to fit the logic behind the 
procedure. Since in the adiabatic approximation a has 



5 



to be a constant, not just a slow variable, this is clearly 
what we have to do. But this point is indeed a problem: 
it is the problem discussed in the next Section. 

Even if this subtlety did not bother us, when the ad- 
justments ( |l7|) in initial conditions and a are made in 
evolution under Hqem and H$s, the first order errors of 
Fig. 4 are reduced, but only to second order; see Figs. 5 
and 6. Although Hss has no real credentials at this 
point, H GEM is supposed to be accurate at 0{e 2 ). So we 
must explain the second reason for its failure. 
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FIG. 5: The errors in the geometric electromagnetism approx- 
imation, with the corrected initial conditions (0). Parame- 
ters are otherwise the same as in Fig. 4; but note the change 
in vertical scales. The thickness of the A8 curve reflects the 
rapid oscillations in the exact motion, which our horizontal 
scale is too coarse to resolve visibly. 
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FIG. 6: The errors in the heuristic approximation of Schmied- 
mayer and Scrinzi, with the corrected initial conditions ([lTj). 
Parameters are otherwise the same as in Figs. 4 and 5. Note 
the vertical scale. The thickness of the A8 curve reflects the 
rapid oscillations in the exact motion, which our horizontal 
scale is too coarse to resolve visibly. 



V: 'GEOMETRIC GRAVITY' 

The second order problem is more profound than the 
first order problem of initial condtions: it is that Hqem 
is incomplete. The scalar potential of GEM is indeed 
a second order correction; but as shown by Littlejohn 
and Weigert, at second order there also appears a geo- 
metric tensor ■potential: the effective Hamiltonian must 
acquire a term of the form Y^,jkPjPk5 1 (%)■ In analogy 
with 'geometric magnetism', we can think of the metric 
of general relativity, and call this second order effect 'ge- 
ometric gravity'. Converted to our original dimensionful 
units, the geometric gravity term derived in Q is 



H, 
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when in this last step we return to our dimcnsionlcss 
variables, and drop terms of higher order in e. This yields 
the full second order adiabatic effective Hamiltonian 



H 



AD2 
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distinct from both H$s and Hqem- 

(We can already notice a point in favour Had2- it is 
numerically equal to the exact H up to second order in 
e, since with the adjusted initial conditions of ( p7| ) we 
find H AD2 (0) = H(0) + C(e 3 ), and both Hamiltonians 
are conserved under the evolutions they each generate.) 

Where does 'geometric gravity' come from? We can 
see this by re-examining the equations of motion for a 
and a in (|i"o|). We can recall a{t) as given in (16), and 
then by inspection we can also find that 



a(t) 



eape 



sin - 



0(e 2 



(20) 



where r(t) = Jdtr . These results then imply that 



3 \/\ — a 2 cos a = s yl 



o-pe 



0(e 2 ). 
(21) 

So in fact the first order correction to the purely fast 
zeroth order term in a z contains the slow term eapg/r. 
Thus at first order in e, it is no longer precisely the spin 
component in the B direction which is the adiabatic in- 
variant. The axis about which the spin precesses rapidly 
tilts slightly in the z-direction, so that the spin compo- 
nent which varies slowly is not exactly sg, but rather 



s ' 



sg 



r 2 I'u 



(22) 
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See Fig. 7. 




FIG. 7: The magnetic particle's spin (two-toned arrow) pre- 
cesses rapidly around an axis (dark arrow) tilted slightly z- 
ward from the ^-direction (grey arrow), which is tangential to 
the magnetic field line (dashed arc). A post-adiabatic approx- 
imation that includes only geometric electromagnetism treats 
the spin component along the grey arrow as an adiabatic in- 
variant. Incorporating 'geometric gravity' as well means us- 
ing a more accurate adiabatic invariant, the spin component 
along the dark arrow. 

This means that the adiabatic approximation is not 
exactly to average over a, but to average over the variable 
conjugate to sgi . If we now re- write the Hamiltonian H 
in terms of cr' = sqi / s and its conjugate angle a', we find 



H 



V r Pa + £ (1 — c ) cos a a' e 

y + ^ w - 

p 2 p 2 g + £ 2 (1 - a' 2 ) cos 2 a' 
2~ + 
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(23) 



Averaging (p3) over a' now produces Hadi-, with <j — > <j'. 
While cr — * er' is really a trivial change, since in the effec- 
tive theory this quantity is simply a constant, this does 
resolve the concern we raised in the previous Section, im- 
mediately after Eqn. (Jl7|). The constant component a of 
o~(t) is the entire slow component of a' (up to corrections 
of third order or higher); in replacing ct'(O) — * tr, we are 
indeed discarding only fast components. 

(An extra complication in this derivation of Hadi is 
that the change of variables from a, cr — > ct',a', with- 
out also changing p r and pa, is not canonical. The non- 
vanishing Poisson brackets [<j',p r ] and [cr' ,9] are purely 
fast quantities, however, and this means that we can ig- 
nore this complication. If we also adjusted the momenta 
to keep our variables exactly canonical, extra terms not 
shown in (|2^) would appear when H was re-written; but 
these extra terms would then vanish in the averaging over 
a'.) 
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FIG. 8: The errors in the correct second order adiabatic ap- 
proximation, evolution under Had2- Parameters are again 
e = 0.01 and the initial conditions on the exact motion are 
the same as in Figs. 4-6, while the initial conditions for the 
Had2 approximation are given in terms of the exact initial 
conditions by Eqns. (^). Note the change in vertical scale 
in comparison with Figs. 6 and 7. On this vertical scale the 
rapid component of the exact motion is quite large, but the 
unchanged horizontal scale is still too coarse to resolve it, 
producing the thick black bands. 



As illustrated in Fig. 8, when e is small, Hadi with 
(|r7| ) is accurate for the slow motion even at Oie 2 ). The 
power of the adiabatic approximation is truly remark- 
able: it ignores the small, rapid 'wobble' in the exact 
evolution, and it tracks the slow motion very accurately 
over many orbits. The contrast with Fig. 4 shows pri- 
marily the effects of correcting the initial conditions ac- 
cording to ( [l7| ) , which involves a first-order change in the 
adiabatic invariant cr from the instantaneous exact <r(0). 
Comparing Figs. 5 and 6, however, shows that even with 
the improved initial conditions, the simpler approxima- 
tions Hqem and Hss are not as accurate as Had2 in 
the illustrated case. We now understand the failings of 
Hqem', and it was not clear in the first place why one 
should expect Hss to be a good effective Hamiltonian. 
But in the case shown in Figs. 1-3, Hss was really ex- 
cellent. Why? 



VI: PREFERRED CANONICAL CO-ORDINATES 

The uncanny effectiveness of Hss, even though it does 
not seem to be the correct adiabatic effective Hamilto- 
nian H ADii is very simply explained: Hss actually is 
the correct adiabatic effective Hamiltonian, in a different 
set of variables. If we perform the canonical change of 
variables r — > r — \s 2 cr in Had2, then to second order 
in e we obtain Hss- Hence the heuristic evolution of 
Schmiedmayer and Scrinzi would disagree with the cor- 
rect adiabatic result only at third or higher order in e, 
if in comparing the two evolutions one also applied the 
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small constant shift in r, and used the more accurate 
initial conditions (|17|). See Fig. 9. 
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FIG. 9: The errors in the heuristic approximation of Schmied- 
mayer and Scrinzi, with the corrected initial conditions (J(LT\), 
and the co-ordinate shift described in the text. Evolution un- 
der Hss proceeds from the shifted initial condition f(0) = 
1.5+ §£ 2 o", and then f (t) — r(t) — \£ 2 o is plotted as Ar. Note 
that the vertical scales are the same as in Fig. 6. 

Using the trivial initial conditions instead of (|l_ 
amounted to an error in a of — cr 2 (0)/r(0), which 
was of order 1% in the case of Figs. 1-3, and not clearly 
discernible in the large-scale plot shown in our first three 
Figures. In the smaller r(0) case of Fig. 4, close ex- 
amination revealed the error. Taking the correct initial 
conditions, but failing to include the co-ordinate shift in 
r, led to the 0(e 2 ) error seen in Fig. 7. Yet this simple 
co-ordinate shift is in a sense the main reason why Hss 
does so much better than Hq em in the strongly process- 
ing case of Figs. 1-3: since the canonical transformation 
relating Hss and Had2 is a constant shift in radius, even 
the most naive calculation with Hss gives the orbital pre- 
cession with full second-order post-adiabatic accuracy. 

Just when everything at last seems clear, though, there 
is one final twist in the classical part of our story: see 
Fig. 10. In some respects the canonical co-ordinates 
yielding Hss as the adiabatic effective Hamiltonian are 
significantly better than those yielding Had2- For one 
thing, Hss itself is more easily solved than Had2, be- 
cause it is effectively just the Kepler Hamiltonian with a 
renormalized pg . 

More importantly, though, Had2 contains a term pro- 
portional to e 2 r -3 , while Hss has no terms more singular 
than the r~ 2 centrifugal barrier. As long as r remains of 
order e°, the more singular term causes no problems; but 
for r < the formally e 2 correction in Had2 is only 
smaller than the zeroth order a jr term by 0(e), and the 
adiabatic hierarchy will begin to break down. Of course, 
since the exact Hamiltonian has nothing worse than r~ 2 , 
one of the things this breakdown implies is that the r~ 3 
term is not really accurate when r is this small. So while 



the r~ 3 term is quite accurate at r > O(e ), its singu- 
larity at small r is spurious - an artifact of an adiabatic 
expansion applied beyond its regime of validity. 

In contrast, no hierarchy breakdown problems occur 
with Hss until r < e\ and even then there is no severe 
disaster, because the zeroth order centrifugal barrier still 
dominates. So the variables of Hss are better than those 
of Had2 in that they greatly shrink the phase space re- 
gion where the adiabatic approach fails, and also reduce 
the severity of the failure. And this advantage can trans- 
late into significantly more accurate results, as shown 
in Fig. 10. The highly eccentric orbit in this case, and 
the relatively large value of e, mean that as the particle 
approaches the wire Had2 suffers significantly in com- 
parison with Hss- 




r cosO 



FIG. 10: The exact trajectory as in Fig. 1 (solid curve), with 
the adiabatic trajectories as given using Hss (long dashes), 
and using Had2 (short dashes). In both approximations 
the initial conditions are given by ( |l7|) . For Hss evolution, 
r(0) = 10 + £ 2 <t/2, and then r(t) -~?a/2 is plotted as the 
radius. This shift, and the improved initial conditions, make 
this corrected implementation of the Schmiedmayer-Scrinzi 
approximation even better than the naive implementation 
shown in Fig. 3: it is barely distinguishable from the exact 
curve. Although the two dashed trajectories should formally 
differ only by order e 3 , the breakdown of the adiabatic hierar- 
chy as the particle approaches the wire has a much more severe 
effect in the Hadi co-ordinates than in the Schmiedmayer- 
Scrinzi co-ordinates. 

This raises an interesting question in adiabatic the- 
ory. When the separation of time scales does not apply 
over the entire phase space, it may be that the nonadia- 
batic region (in which there is no time scale separation) 
can be patched into the adiabatic region (in which there 
is good separation of time scales) more smoothly with 
some canonical coordinates than with others. Is there a 
systematic way of identifying the better co-ordinates, in 
general? In this paper, we can do no more than raise this 
question as a subject for future investigation. 

In classical mechanics this issue is not a problem if one 
simply avoids the nonadiabatic region, such as by consid- 
ering orbits that maintain r >> s. In quantum mechan- 
ics, though, the problem is in principle much more severe. 
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The energy spectrum of Hss is well behaved, and can 
even be obtained exactly (though Hss itself is of course 
only a second order approximation); but the spectrum 
of Had2 with positive a is continuously unbounded from 
below. One can of course recognize the spuriousness of 
these negative energy states, and try to discard them in 
various ways. (Semiclassical calculations, for instance, 
will essentially reduce the advantages of Hss over Hadi 
to only those that we have noted classically.) If numerical 
methods are required, however, as they may well be for 
more complicated problems than the one we have been 
considering, a spurious singularity introduced by adia- 
batic elimination might be quite a serious obstacle. 



VII: THE QUANTUM VECTOR KEPLER 
PROBLEM 



the scalar Schrodinger equation 
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(27) 



which is clearly the quantized version of Had2- 

If we at last define • v /r$(r) = \frf(r — e 2 a/2) to ef- 
fect the change of variables from Had2 to Hss without 
problems from the r~ 1 d r term, we obtain the quantized 
version of Hss- 



E f = ~^F> 



2J 2 



dr 2 



19/ 

r dr 

- a -f- 



l + i[ £ (£+f)-3e 2 a 2 l 



2r 2 



/ 



(28) 



That such problems can indeed arise quantum mechan- 
ically is confirmed by quantizing our vector Kepler prob- 
lem. Keeping track of what our rescalings have done to 
the commutation relations, we find the time-independent 
Schrodinger equation corresponding to (|^): 



1 / h z 

= r 

2 V J 



dr 2 



ld_ 

r dr 



(1 -ea z f 



r 



(24) 

where ^ = ty(r) is spinor- valued, a x z are the spin-s 
generalizations of the Pauli spin matrices divided by the 
total spin s, and E is the dimensionless energy eigenvalue. 
To obtain the equation in this form we have already used 
a (9-dependent spin basis such that s x is in the local 9- 
direction, and set pe — ► 1. 

Remembering what we learned classically about the 
tilting of the adiabatic invariant, we can make a short- 
cut to the correct approximation at second order in e by 
introducing 



(25) 



Since h/J is certainly of order e, the terms that arise 
from the commutator between this r-dependent rotation 
and the kinetic term are of 0(e 3 ); but we should still 
keep the term linear in d ri since {%/ J)d r ~ p r may be of 
order £ . So up to second order, the result is 



Eil) = 



J 2 

-is 



dr 2 
h 2 g y d 
J 2 ' 



ld_ 

r dr 



1 



2 g r 



r 



1 



e 

2T 2 



xjj . (26) 



The quantum analogue of the classical adiabatic ap- 
proximation for problems like this one is the Born- 
Oppenheimer approximation. In ( p6| ) we now implement 
the Born-Oppenheimer approximation by considering an 
eigenstate of cr x , with eigenvalue a. We take the expec- 
tation value of the Hamiltonian in this state, and solve 



In principle then we must impose boundary conditions at 
r = £ 2 er/2 instead of at r = 0, but since at this small r we 
will have / ~ r J / n , we may safely ignore this point. The 
eigenvalue equation for / may then be solved exactly Jl0[, 
giving us an adiabatic approximation to order e 2 for the 
bound state energy eigenvalues of the exact problem: 



a J 



1 



2fr 2 (n 



mj) 



(29) 



where n £ [1,2,3,...], and the 'adiabatically corrected' 
angular momentum quantum number is 



1 



mj =-y/j 2 + [s{s + h)- 3s 2 a 2 ]/2 



(30) 



Note that ( p9| ) is also obtained, to order e , by WKB 
quantization of (|l9|). 

For the special case s = h/2, (E9|) is not only accurate 
to C(e 2 ), but is actually exact [jllj. This was precisely 
what motivated Schmiedmayer and Scrinzi to consider 
their effective Hamiltonian for general s. For other val- 
ues of s and a range of values of J, numerical results in 
Ref. H show that ( p9| ) is actually accurate to 0{e 4 ). The 
success of the heuristic formula in ||] appeared to be the 
success of a mysterious rival to standard adiabatic the- 
ory; but the mysterious rival has now been revealed as 
standard adiabatic theory in a thin (but evidently useful) 
disguise. 

Comparison with most other treatments of the quan- 
tum vector Kepler problem is not relevant, because they 
examine predominantly cases of small J, where the adi- 
abatic method does not perform well. Ref. |l2| provides 
numerical results in only one case where e is really small, 
namely the exactly solvable case s — h/2 and J = 117i/2. 
The reported results of the adiabatic calculation in E^] 
show errors of order e 4 in this case. As far we under- 
stand, it is only fortuitous that our results are actually 
better than this; but perhaps some advantage is due to 
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the fact that the adiabatic technique used in |12[ is un- 
controlled, whereas ours is a systematic expansion in the 
small parameter e. The effective potential (for the bound 
channel) used in jl^ ] is, in our units, 

1 + e 2 , 

U- = —V " Vr- 2 +e 2 r-4 , (31) 

which is the exact lower eigenvalue of the the 2-by-2 po- 
tential matrix of ( |24| ) for s = h/2. Expanding to C(e 2 ) 
produces the potential Had2, with its problematic r~ 3 
behaviour. As we have seen, the canonical transforma- 
tion to Hss effectively removes this spurious singularity. 
Since ( pl| ) is also well-behaved as r — > 0, and has the 
correct behaviour to 0(e 2 ) at large r, it is not surprising 
that it also gives results that are good at 0(e 2 ). Ob- 
taining still better results, however, requires more than 
simply expanding {/_ to higher order in e. 

VIII: DISCUSSION 



that the adiabatic effective Hamiltonians may have sin- 
gularities like the r -3 potential of our Had2, which are 
unphysical artifacts of applying the adiabatic formalism 
beyond its domain of applicability. 

In difficult cases of this kind, whether quantum or 
classical, a kind of connection formula would seem to 
be needed, from which one could compute the time and 
phase space location at which the system would exit the 
region of adiabatic breakdown. In the problem studied 
in this paper, however, we have seen that a co-ordinate 
change can eliminate a spurious singularity, producing 
an effective Hamiltonian which does not seem to need a 
connection formula. It would clearly be desirable to de- 
termine whether this is simply a fortunate coincidence 
in this one case, or whether there may exist a general 
method of identifying the best co-ordinates, and mini- 
mizing the damage done by what we have called crisis 
regions. 



Adiabatic methods are a remarkably powerful analyt- 
ical tool. Their power often does seem to be remark- 
able, in the sense that they frequently deliver accuracy 
greater than one can naively expect. The robustness of 
adiabatic approximations has perhaps not yet been fully 
understood. For example, an interesting issue that we 
have raised in this paper concerns the possibility that 
the hierarchy of time scales, upon which adiabatic meth- 
ods depend, may break down within a localized region of 
phase space. A classical system may then evolve outside 
this region for most of its history, so that its behaviour 
can be closely approximated adiabatically, but in a suc- 
cession of intervals of 'crisis', it may pass briefly through 
the region where adiabatic methods fail to apply. During 
the interval of crisis the fast degrees of freedom do not de- 
couple from the (nominally) slow ones, and so the precise 
state of the fast variables might affect the slow variables 
in ways that could persist after the crisis. Thus the fast 
degrees of freedom, which within the adiabatic regime are 
effectively hidden variables, could 'emerge' during crises. 

While this kind of emergence of hidden variables has 
(presumably) nothing to do with quantum measurement, 
the classical phenomenon will have analogues in quantum 
mechanics. The effective Hamiltonians given by standard 
adiabatic techniques may be unreliable in some regions 
of Hilbert space. The result may in principle be spuri- 
ous simplicity, but it may also be spurious complexity, in 
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